Entanglement and errors in the control of spins by optical coupling 



in 

O 

o 



5h 

43 



E 

C 

o 
o 



> ■ 

-xt- : 
^t- . 

m ■ 
O ■ 
in : 
o . 

-I— > ' 

a ■ 

E 

T3 ■ 

c : 
o . 
o ■ 



X 



QD (resonant) 


Light 


QD 

• (off resonance) 








V- 








Impurity spins 



G. F. Quinteiro and C. Piermarocchi 
Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824 

(Dated: February 2, 2008) 

We analyze the optical quantum control of impurity spins in proximity to a quantum dot. A 
laser pulse creates an exciton in the dot and controls the spins by indirect coupling. We show 
how to determine the control parameters using as an illustration the production of maximal spin 
entanglement. We consider errors in the quantum control due to the exciton radiative recombination. 
The control errors in the adiabatic and nonadiabatic case are compared to the threshold needed for 
scalable quantum computing. 



I. INTRODUCTION 

Impurity spins embedded in semiconductors are cur- 
rently under investigation for quantum computing imple- 
mentations. Recently, optical techniques have been pro- 
posed to control the spin-spin coupling and realize two- 
qubit quantum gatesi 1 ^ 3 ^ 5 ^ 7 The optical method sug- 
gests the possibility of an ultrafast control of the qubits. 
The flexibility in the control that can be obtained by 
pulse shaping 8 and the absence of noisy contacts repre- 
sent additional advantages. On the experimental side, 
ensemble optical measurements have demonstrated the 
production of spin entanglement for impurities embed- 
ded in a semiconductor host£ More recently, the mea- 
surement of the quantum state of a single impurity spin 
obtained by coupling it to a single exciton in a quan- 
tum dot (QD) has been experimentally carried out^S In 
this paper we study theoretically the control of impurity 
spin states when the interaction among them is controlled 
by optically-generated excitons in a QD. We analyze the 
control errors due to the radiative recombination of the 
exciton that mediates the interaction between the spins. 
Moreover, we illustrate how the control parameters can 
be obtained directly from simple analytical expressions. 
The method is applied to design the control parameters 
in the production of maximal spin entanglement. 



II. SYSTEM 

The physical system consists of two impurity spins 
placed close or inside a QD in such a way that there 
is not a direct interaction between them. A schematic 
picture is given in Fig. ^ By coding the qubit in more 
than one spin efficient schemes for fault-toleranljii and 
exchange-onljii^ quantum computation can be naturally 
applied to this setup. Dots of different size provide the 
frequency selectivity to address specific spin pairs and 
realize two-qubit readouts. The model we use contains 
few parameters describing the exciton-light and exciton- 
impurity coupling and can be applied to different physical 
systems. For instance, it can be used for excitons local- 
ized by monolayer fluctuation in III-V and II- VI quantum 
wells and interacting with a finite number of localized 
impurities as in Rcf. 9. III-V or Si/Ge self-assembled 



FIG. 1: Scheme the system: localized spins located near self 
assembled QDs are coupled by an exciton created by a laser 
pulse. Dots of different size provide selective control and read- 
out. 



QDs can also be used as shown in Fig. ^ For typical 
semiconductor systems we can restrict the analysis to 
heavy-hole excitons due to the splitting between heavy- 
hole and light-hole bands in the dot. The heavy-hole 
exciton spans a four dimensional space consisting of two 
optically-active and two dark states. We treat the inter- 
action between the electromagnetic field and the excitons 
semiclassically, and we consider spin states that inter- 
act only with the photoexcited electron in the dot. This 
is the case for instance of donor impurity spins in typi- 
cal semiconductors because the electron-hole exchange is 
much smaller than the electron-electron exchange. No- 
tice that by using circularly polarized light the exciton 
induces, besides the spin-spin coupling, also a local effec- 
tive magnetic field on the spins. 1-13 This effective mag- 
netic field can be controlled by the laser polarization and 
disappears for linearly polarized light. We will consider 
below the case of circularly polarized light. 

The exciton-spin part of the Hamiltonian can be writ- 
ten as 



H = eN - 2V(J 2 - S 2 - I?) 



(1) 



where S is the total spin of the two impurities (h = 1) and 
e is the energy of the exciton in the dot. The operators N 



and L are defined as N = bib 
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FIG. 2: Energy level diagram and optical selection rules for 
(j+ polarized light. In n = 1, the total splitting is j = 6V. 

and hole spin +3/2, while &j creates a dark state exci- 
ton with electron spin +1/2 and hole spin +3/2. The 
total angular momentum Ji = Li + Si, and V is the ex- 
change interaction between the impurity spins and the 
photoexcited electron in the dot. The strength and the 
sign of V depend on the system. For instance, this cou- 
pling is expected to be ferromagnetic for electrons in the 
dot interacting with localized rare-earth magnetic impu- 
ritites, while it is antiferromagnetic for a dot mediating 
the interaction between shallow donorsi^ Without loss 
of generality, we will assume V > below. The coupling 
of the excitons in the QD and the external laser field is 
given by 

H (LX) = l) e -i,i 6 t + hc (2) 

2 1 

where f2(t) is the time-dependent Rabi energy associated 
with the optical pulses, and ix> is the energy of the laser. 
We consider only anti-clockwise polarization (o + ) which 
generates excitons with electron and hole spin states 
— 1/2 and +3/2, respectively. Excitons with hole spin 
—3/2 are not included in the model since they are not 
excited by a + light and the impurity spins can only flip 
the spin of the photoexcited electron. 

A scheme of the relevant energy levels is given in Fig. 
12 In the ground state n = we have the singlet and the 
triplet states corresponding to the two non-interacting 
impurities. In the excited state n = 1, the electron in 
the dot splits the triplet states in a quadruplet J = 3/2 
and a doublet J = 1/2. The total Hilbert space is thus 
spanned by a total of 12 states. The arrows in the scheme 
identify the selection rules for cr+ optical transitions. The 
transitions have different oscillator strengths, which are 
calculated using the Clebsch-Gordan coefficients. Notice 
that the light does not connect directly states with differ- 
ent spin S. The structure of the energy levels provides a 
natural readout scheme for the coded logical qubit |0l), 
1 1 l) in the exchange-only schemed An optical setup sim- 
ilar to the one for single spin readout^ could be used: a 
single peak at e corresponds to |0l) while two peaks sep- 
arated by 6V correspond to the logical state |1l)- 



III. QUANTUM CONTROL 

In order to illustrate how to design the optical con- 
trol we consider the production of maximal spin entan- 
glement. We choose the the initial state | TO ® |0) as 
the tensor product of a linear superposition of impu- 
rity states {IS 1 ), |T )}, and the exciton |0) representing 
an empty QD. We consider separately the case of infinite 
and finite 7 , i.e. spontaneous radiative recombination 
lifetime for the exciton in the dot. In the first case we 
determine analitically the control parameters that pro- 
vide maximal spin entanglement. In the second case, we 
solve numerically the master equation for the full system 
in Fig. [3 This will allow us to analyze errors due both 
to the radiative recombination and to the finite proba- 
bility of remaining with an exciton in the dot at the end 
of a pulse. The latter is an error similar to a double 
occupancy error in the case of spins controlled by gate 
voltagesi*. Ideally, the QD must be empty at the end 
of each optical pulse, and this can be achieved by an 
adiabatic evolution, or by a nonadiabatic evolution plus 
additional conditions in the pulse areai^. 



A. Infinite radiative lifetime 

We first analyze the ideal case of a nonadiabatic evo- 
lution at 7 = 0. We call nonadiabatic the evolution that 
follows from a laser resonant with at least one transition 
between the n = and n = 1 subspaces in Fig. |3 This 
implies that there is a substantial exchange of energy 
between the electromagnetic field and the dot, which in 
turn results in a significant population inversion during 
the pulse. Using a numerical simulation we illustrate in 
Fig.|3]the evolution of the \S) and |T ) populations under 
a Gaussian pulse giving a Rabi energy of the form, 

m = ^e-^f . (3) 

The pulse is resonant with the bare exciton energy, which 
in the scheme of Fig. [21 corresponds to a resonant transi- 
tion for the singlet state. In order to have no excitonic 
population at the end of the pulse, we need the pulse area 
for the resonant excitation to be multiple of 27r, therefore 
f2 and r are chosen so that the pulse area is J7r = 8tt. No- 
tice that the population of the ground state singlet \S) 
is completely depleted during the pulse but at the end 
comes back to the original population (0.5). In contrast, 
the triplet (|To)) population follows an adiabatic evolu- 
tion due to the exchange interaction affecting the optical 
resonance. In Fig. [3] (inset) we show the real and imag- 
inary part of the coherence (iS|j(?|To). In order to create 
the maximally entangled state we need a ±7r/2 phase in 
this matrix element and the chosen optical pulse achieves 
this goal. This relative phase transforms, for example, 
the state | U) into 2- 1 / 2 (|S') + i\T )) oc | U) + *| IT)- 
For a given value of the exchange coupling V and pulse 
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width t, the maximum intensity of the field in Eq. © 
is found from the roots of the equation 



4>t(SI,V,t) ±n/2 = 



(4) 



where (f>r = f_ A^ (t) dt is the dynamic phase that the 
state \Tq) picks up following the adiabatic evolution. No- 
tice that since the pulse is a multiple of 2ir the singlet 
will only pick up a trivial phase (±1). A^ is the eigen- 
value satisfying A^(±oo) = for a 3-level Hamiltonian 
representing the triplet states, 



H T {t) 
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(5) 



The optical detuning 5 — e— lo and j = 6V is the splitting 
in the excited state between J = 3/2 and J = 1/2 states. 
If we assume that the three eigenvalues of the matrix 
in Eq. [S] do not cross during the pulse evolution, the 
expression for A^(t) can be written as 



3 V 3 



where 



with 



6(t) = 2kn- 



arccos 



r(t) 



(6) 



(7) 



q(t) = 

r(t) = 4(4j- 



+ 12 j5 + 365 2 + 27Q 2 {t) 

38){5j + 35) {j + 65) - 81(2j + 35)n 2 (t) 



If the exciton impurity coupling is ferromagnetic {j > 0), 
we have to take in Eq. (J7| k = 1 for <5 > j/3, k = 2 
for -2j/3 < 5 < j/3, and k = 3 for 5 < -2j/3. In 
contrast, for j < 0, we have to take in Eq. J7J) k = 1 for 
5 > -2 j/3, k = 2 for l/3j < 5 < -2j/3, and k = 3 
for 5 < j/3. The analytical expression in Eq. JBJ) allows 
us to determine exactly the control parameters from the 
roots of Eq. (gj . 

In the 7 = adiabatic regime the laser pulse is tuned 
away from the optical resonances between the n = and 
ri = 1 levels. An example of a simulation of an adiabatic 
control is shown in Fig. The laser is tuned 2 meV 
below the bare excitonic energy and 1 meV below the 
triplet resonances corresponding to J = 3/2 in Fig. [21 
We plot in Fig. 0] the same quantities of Fig. |3J Notice 
that in this case the pulse area can be arbitrary, provided 
the adiabaticity is preserved. The change of phase in 
the coherence (S*|/9|To) is now obtained with a smooth 
transition. The control parameters in this adiabatic case 
are determined by the roots of 



6r(fi,V,T)-&j(fi,V,T)±7r/2 = 0. 



(8) 
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FIG. 3: (Color online) Nonadiabatic control, 7 = 0. Evo- 
lution of the \S) (solid blue line) and |Tb) (dashed red line) 
populations under a Gaussian pulse of area 87r. The temporal 
width of the pulse r is 7.02 ps and the ratio fi/V is 0.6697. 
(Inset) Real (dashed green line) and imaginary (solid red line) 
part of the coherence (S\p\To). 



0.25 — 



-t 



0.25 




25 



50 75 
time (ps) 



100 



FIG. 4: (Color online) Adiabatic control, 7 = 0. The laser is 
tuned in 2 meV below the bare excitonic energy. Evolution of 
the \S) (solid blue line) and |To) (dashed red line) populations 
under a Gaussian pulse. The temporal width of the pulse r 
is 10.2 ps and the ratio fl/V is 1.24. (Inset) Real (dashed 
green line) and imaginary (solid red line) part of the coherence 
(S\p\T ). 



In contrast to the case of Fig. the singlet now picks up 
a nontrivial dynamic phase (j>s — J-oa Ag(£) dt where Ag 
is the eigenvalue of the singlet Hamiltonian 



Hs(t) = \ 



Sl(t) 
fl(t) 25 



(9) 



with the property Ag(±oo) = 0. As for Ht this has a 
simple analytical form \%{t) = ^±\^/5 2 + fl 2 (t) , (+ for 
5 < and — for 5 > ) which can be used to determine 
the control parameters from the roots of Eq. JHJ . 

B. Finite radiative lifetime 

In order to determine how this control scheme is af- 
fected by the finite lifetime of the exciton in the dot we 
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introduce a finite value for 7, and solve the master equa- 
tion p = — i[H, p]+£[p] using the values of the control pa- 
rameters corresponding to the evolution of Figs.|3and21 
C[p] is the Liouvillian superoperator that can be written 
as £[p] = L^pL + ^{L^L,p} where = ^/j accounts 
for the spontaneous radiative recombination of the exci- 
ton in the dot. Once p is obtained, a 4 x 4 reduced density 
matrix for the impurity spins ps is computed by tracing 
out the exciton degrees of freedom. The entanglement in 
the Bell state is mostly sensitive to decoherence processes 
and its analysis provides a good test for the scheme. We 
quantify the error on the reduced density matrix ps using 
two different methods, the Purity and the Peres criterion 
of separability^ According to the Peres criterion a state 
is entangled iff E m i n < 0, where E m i n is the minimum 
eigenvalue of a matrix constructed by transposing the 
non-diagonal 2x2 blocks of ps- A maximally entangled 
state has a E m i n = —1/2. The deviation from that value 
gives a measure of the effect of the radiative recombina- 
tion on the entanglement and we quantify the entangle- 
ment error as AE — E m i n + 1/2. The purity of ps is 
a different parameter that characterizes the error in the 
spins states due to their entanglement with the exciton 
in the dot. We quantify this error as AP = Trp 2 s — 1. In 
principle there are errors that can disentangle the spin 
states without a change in the purity, for instance by 
affecting the phase picked up in Eq. JSJ. Therefore, in 
principle the errors induced by 7 affect independently 
AE and AP. 

We compare in Fig. [3] the errors in the Entanglement 
AE (upper) and in the Purity AP as a function of 7 for 
both the adiabatic and nonadiabatic evolution. Both er- 
rors increase linearly at small 7. However, the errors in 
the adiabatic case are always smaller than in the nona- 
diabatic case in the range of parameters we have inves- 
tigated. We remark that, due to the incommensurabil- 
ity of the eigenvalues of Hs and Ht, there are not spe- 
cial conditions that would give perfect entanglement with 
square pulses as in the case of a direct spin-spin cou- 
pling^ An important figure of merit for the application 
of this quantum control technique to quantum computa- 
tion is provided by the error per gate parameter. This 
has to be below a threshold value in order to make scal- 
able quantum computing possible. The estimate for such 
a threshold depends on assumptions on the error model 
and device capabilities but the 10 -4 value — is usually 
used as a benchmark in typical experimental implemen- 
tations. The error in the entanglement production gives 
an estimation of the error per gate since the quantum 
operation done corresponds to a V SWAP modulo some 
single qubit operations. We see in Fig. [3] that the 
threshold can be achieved for 7 smaller than IpeV. Self 
assembled QDs have typically a ground state exciton life- 
time of the order of one or more nanoseconds and would 



reasonably be in this region of parameters. 



Log(AE) 




FIG. 5: (Color online) Log-Log plot of the deviation from 
maximal entanglement (upper panel) and maximal purity 
(lower panel) as a function of the radiative recombination 7 
(in meV). Solid blue line: adiabatic evolution. Dashed red 
line: nonadiabatic evolution. 



IV. CONCLUSIONS 

In conclusion, we have analyzed the entanglement pro- 
duction between two spin-impurities induced by an exci- 
ton in a neighboring quantum dot. In the case of 7 = 0, 
the parameters for the quantum control can be analyti- 
cally determined from the roots of simple integral equa- 
tions. We showed that the finite lifetime 7 -1 of the exci- 
ton in the dot can affect the purity of the spin states and 
introduces errors in the entanglement production. In ad- 
dition we found that such errors increase linearly with 7 
and can be kept below the 10 -4 threshold for error cor- 
rection if parameters typical of self assembled QDs are 
used in the simulation. 



Acknowledgments 

This work was supported by the NSF Grant DMR- 
0312491. We thank Prof. T. A. Kaplan for discussions. 



1 C. Piermarocchi and G. F. Quinteiro, Phys. Rev. B. 70, 2 C. Piermarocchi, P. Chen, L. J. Sham, and D. G. Steel, 
235210 (2004). 



5 



Phys. Rev. Lett. 89, 167402 (2002). 

3 R. Rodriquez, A. J. Fisher, P. T. Greenland, and A. M. 
Stoneham, J. Phys. Condens. Matter 16, 2757 (2004). 

4 A. M. Stoneham, A. J. Fisher, and P. T. Greenland, J. 
Phys. Condens. Matter 15, L447 (2003). 

5 G. Ramon, Y. Lyanda-Geller, T. L. Reinecke and L. J. 
Sham, |ajXrvxond -mat/0412003 (2004) 

6 A. Nazir, B. W. Lovett, S. D. Barrett, T. P. Spiller, and 
G. A. D. Briggs, Phys. Rev. Lett. 93, 150502 (2004). 

7 E. Pazy, E. Biolatti, T. Calarco, I. D'Amico, P. Zanardi, 
F. Rossi, and P. Zoller, Europhys. Lett. 62, 175 (2003). 

8 P. Chen, C. Piermarocchi and L. J. Sham, Phys. Rev. Lett. 
87, 067401 (2001). 

9 J. M. Bao, A. V. Bragas, J. K. Furdyna, and R. Merlin, 
Nature Mat. 2, 175 (2003); J. M. Bao, A. V. Bragas, J. K. 
Furdyna, and R. Merlin, Solid State Commun. 127, 771 
(2003); J. M. Bao, A. V. Bragas, J. K. Furdyna, and R. 



Merlin, Phys. Rev. B 71, 045314 (2005). 

10 L. Besombes, Y. Leger, L. Maingault, D. Ferrand, H. Ma- 
riette, and J. Cibert, Phys. Rev. Lett. 93, 207403 (2004). 

11 D. Bacon, J. Kempe, D. A. Lidar, and K. B. Whaley, Phys. 
Rev. Lett. 85, 1758 (2000). 

12 D. P. Di Vincenzo, D. Bacon, J. Kempe, G. Burkard, and 
K. B. Whaley, Nature 408, 339 (2000). 

13 Combescot and O. Betbeder-Matibet, Solid State Com- 
mun. 132, 129 (2004). 

14 X. Hu and S. Das Sarma, Phys. Rev. A 61, 062301 (2000). 

15 T. A. Kaplan and C. Piermarocchi, Phys. Rev. B 70, 
161311(R) (2004). 

16 A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

17 D. Gottesman, Stabilizer Codes and Quantum Error Cor- 
rection PhD thesis, Calif. Inst. Tech. Pasadena, California 
(1997). 



